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k tiaplt to dofcribo fast plaMt haatlog b;f aaeaalotti and iMrtial 

raaiatlvlty affacta la praaantad. A aaall fraction of tha plaaaa 
eontaina atrong eurranta chat cun par alia 1 to tha aagnatlc flald and ara 
dcivan by an axpoeantiaclng alaetclc flald. 'Sta woaaloua eharaetar of tha 
eurrani ^laalpatlon la cmaad by tha axcltatlon of alactcoataClc lon*eyclotron 
and/or lon*Mouacle vavaa. Tha peaalbla cola of raalatlvlty dua to gaooatrlcal 
affacta ('loactlal raalatlvlty') la alwo emaldarad. Through tha uaa of a 
■ar^aal atd>lllCy analyala, aquatlona for tha avaraga alactron and Ion taa- 
paraturaa ara darlvad and autBarlc'aily solved. No loas oachanlsas have baan 
CaiMn Into accent. Tha avolutlon of tha plaaaa la dascrlbad as a path In tha 
drift valoclcy dlagraa. In which tha drift velocity v^/v^ is plotted as a 
function of the electron to ion taaparatura ratio T^/Tj^. 

For currant l^ars with dlMnslons chat ara large coaparad to tha affective 
electron aaan free path, Inertial resistivity Is not Important *and find 
that: 

(1) Beating dua to classical resistivity cannot make T^/T|, >3.1 before 
Instability sets In. This Is always tha lon-cyelotron Inst^lllty. 

(2) father or not a situation with T,/T^ » 1 will occur depends critically 
on tha saturation level of ion-cyclotron waves. 

(3) Ion-acoustic wave heating produces a limiting X^/T^ of 6.4. 

If Inertial raslativlty is Important much hi^ar values of Xg/Tj^ can be 
atcalnad, because inertial resistivity priaarily affects electrons. 

Tha assumption that hard X-rays, aaittad during tha impulsive phase of 
a solar flare, are thermal in origin, requires a hot thermal plasma, 

,Tg 5el0^ X. Our results indicate proaislng possibilities for the production 
of such a hoc plasma. 

Kay words: plaaaa heating - current dissipation - marginal stability - 


solar flare. 



mraoDocTion 


la a plataa. eoavarslon of aaonotic anorgy into haat can taka plaea by 
currant diaaipation. Oftan tha haatlng proeaada aloM^ao that Couloab ralaxatlon 
batwaan alaetrona and Iona kaapa tha taiqparaturaa of both apaclaa aqual. aaomvi 
if tha anargy input oecura auffieiantly faat* and eaaaa to ba aqual, 
alnea tha anacgy diaaipation OMchaniaa favoura anargication of aithar alaetrona 
or iona, and this is «diat wish to atu^ in our praaant work. 

Na hava two apaeifle intacastat 

(a) Xn tha plawaa-astrophyaica litaratura ona oftan asausMa T^>> on tha 
ganaral grounds that oost of tlM dlasipatad anargy goaa into alaetrona 
and tlM alactron-ion ralaxatlon tlaa is ralativaly long (a.g. Kaplan and 

.Tsytovieh, 1973). Na want to aaa if this assuB^tion can ba justiflad and i: 

so, «diat valuas for T /T. can ba attalnad? 

• * ( 1 ) 

(b) k particular axai^la of rapid plasaa haating is tha solar flara. Xt 
has baan found that tha hard X-ray aaission during tha ii^ulsiva phaaa 
of a solar flara can ba attractivaly axplainad as tharaal braasatrahlung 
froa a vary hot plaaos} ^ S 10 X and T^>> (Brown at al. 1979; 
telth and Lilliaquiat, 1979) . Tha authors qvMtad taka tha haating of 
tha plasaa for gran tad and study tha subaaquant avolutlon of tha hot 
raglon aituatad at tha top of a loop, loading to tha formation of a 
conduction front. Wa wish to concantrata on tha aachanisas that ara 
rasponaibla for tha rapid plasaa haating. (Xssarvations indieata that it 
oust occur in a faw aaeonda. Tha paraaatars of tha nuaarical casa atudiaa 
that wa shall praaant ara all tailorad to tha flara haating problaa. 

Rapid eonvarsion of aagnatlc anargy into tharaal anargy* must taka placa in 
local irad ragions, whara tha currant is vary concantratad. Tha raason for 
this is that a larga voluaw with tha aaiM currant dansity avarywhara would 
giva rlaa to unraasonably larga aagnatlc fialds at tha boundary. Also, bacauia 
parallal currants attract ona anothar, thay tand to clus^ togathar (taaring 
instability; Spiear,l9tO; $ehnaek and Killaan, 1978). Wa shall rafar to thasa 
ragions of hi^ currant dansity as 'currant layars'. Thay coincida with tha 
surroundings of nautral linas of tha (fast) tanring aoda configuration and 


(1) Bara and in what follows wa undars^and haating to maan 'bulk anargizaticn' 
Blactron and ion valocity distributions ara assumad to ba roughly isotropic 
and Maxwallian. 
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tlMir largt te«I« apael«l dtatrtbution ew ba ^Ita ehMtlc (Splear* I960). 

Xa th% ewrrane layars eba dlaaipacion rata la ^MraULaad by dia aagaltuda 
of tiia currant danaity tha raaiativlty of tha pXaana. Tha lattar ^aatity 
dapanda on whathar or not tha otrraat la unatabla to tha ganaratlon of vavaa 
(aaa for axaapla Papadopouloa, 1977) and on tha laoaatry of tiM currant layara 
(Inartlal raaiativlty). Va ^all conaldar alactroatatlc Ion-cyclotron vavaa, 
vhl^ haat prlaarily 1 mm, nd lon-aeouatlc vavaa, vhlch haat priaMrlXy 
alactrona. In tha praaanca of thaaa vavaa va calculata tha raaiativlty by 
applying tha narglnal atablllty analyala concapt (MantMlnar, 1977; Manhalaar 
and Borla, 1977). In our analyala va ahall darlva md twoMrlcally aolva tvo 
taaparatura aquationa, ona for tha alactrona and ona for tha Iona. Spaelal 
attantlon vlll ba paid to tha ratio baeauaa thla quMitlty dataralnaa 

vhlch vavatypa haa tha lovar '^hraahold (Klndal and Kannal, 1971). 

A difficult quaatlon la hov tha locally high dlaaipacion rataa around tha 
nautral llnaa Influanca cIm avolution of tha taarlng inatablllty. .. draatle 
alopllflcaClon la naeaaaacy. Va ahall daeoupla tiia haatlng fron tha avolution 
of tha taarlng soda In that va lac tha alactrlc field In tha currant Inyaca 
grov vich a conatant grovch rata . 

XI. MODEL 

Va conaldar a fully lonlaad hydrogan plaama of voluoa V, containing a 
fllaaantary currant dlatrlbuclen. Iha currant layara comprlaa a fraction c 
of tha voluaa V «id In tlMaa currant layara Chara la a atrong and growing 
alactrlc flald In tha fluid raat fraaM, virtually parallal to tha local nag- 
natlc flald. In Uia ranalning fraction 1-c tha alactrlc field in tha fluid 
fraaa la aaauaad to ba affectively aero. Va define V aooMviiat looaaly aa the 
envelope of tha currant Iqyara where rapid dlaaipacion takas place. In our 
dlscusalon the actual aiaa of V ia Isauicarial. Tha dissipation occurs only 
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ia tba fraction c of V and c will ba a small nuobar. Though In ganaral a 
function of tlaa, va shall traat c as a constant. Typically «a shall taka 
c ■ O.Olt without such Justification. 

If tha haat spreads sufficiently fast around the current layers, It Is 
reasonable to Introduce one average electron temperature T^ and one average 
Ion temperature Tj^ In V. Parallel and perpendicular temperatures are taken 
e^al, T^j^ * ^e| ^1<>> " ^ln* assumptions about Tg and T^ require 

some ftirther consideration, for which we refer to Appendix Z. We shall Ignore 
all energy loss mechanlsma from the volume V. These include radiation losses, 
expansion losses and conduction losses. Though Ze Che case of the solar flare 
the latter two will certainly be Important uhen the plasma Is very hot, 
we deliberately Ignore them In this paper, because ve wish to concentrate on 
the heating mechanisms. We plan to Include these loss mechanisms In a later 
study. 

2.1. Dissipated energy end its distribution over electrons and ibcs 

Tha frictional force that determines the resistivity depends very much on 
the state of the plasma. If the resistivity Is caused by Coulomb collisions. 

It depends on the electron temperature only. However when the plasma Is In a 
turbulent state, the electrons carrying the current will also Interact with 
the electric field fluctuations In the waves, which changes the resistivity 
(and othar transport coefficients) of the plasm* in a way that depends upon the 
type of waves that grow. As waves that affect the resistivity we shall consider 

electrostatic lon»cyclotron (IC) and lon-acoustlc (IA) waves (c.f., solar 
applications by Rosner at. al. (1978) and Rinata (1979). 

Another factor of Importance for the resistivity will be the geometry of 
the current layers. Fast electrons crossing smal'. current layers will remain 
under the Influence of the electric field only f . / a short period of time and 
therefore their contribution to the current den Is small. This is due to 
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th« InTClt of ch« eUecronc. The corresponding resistivity is called inertial 
resistivity (Speiser, 1970). 

Inside the current layers several independent scattering processes i occur at 
the same time. Thus, the total dissipation race is related to an effective 

resistivity n^j. »* follows: 

*wt ■ V J • "•''d (1) 


Vq is the drift velocity corresponding to Che current density J. The electric 
field E and j are related through Ohm's law: 


® - V ^ 

where the total resistivity is the sum of individual contributionj: 

^EP ■ J 


< 2 ) 


(3) 


Four reaistivicies are distinguished in this paper (see also Appendix XI5) : 

(a) classical Coulomb resistivity ^CL 

(b) resistivity due to geometrical effects: inertial resistivity 

(c) resistivity due to ion-cyclotron waves: 

(d) resistivity due to ion-acoustic waves: 

Because we are Interested in the behavior of T^ and Tj^, we want to know 

how the dissipated energy is distributed over electrons and ions. Therefore 

2 

we introduce the quantity Xj, ss the fraction of the dissipated energy J 
that goes into the electrons by scattering process i. Classical and inertial res- 
istivities energize only electrons (to order (m/M)^): 


X 


CL 


1 


'IN 


(4a) 
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Cons«>iu«nCly, without w«v«s, th« ions «r« h«mt«d only a« a raault of classical 
aquillbratlon, i.a., by tha' transfer of anargy from alactrons to Ions by 
alactron-lon collisions. Howavar, in tha prasanca of wavas, tha ions can also 
gain anargy from tha wavas. In fact, what happens is that tha waves act as 
an IntaroMdiary, enabling a transfer of anargy from alactrons to ions. Thus, 
despite tha fact that the electric field does work only on the electrons, 

(1 ■ "IC" or "lA") is found to be less than 1. An expression for Xi 1* 
(Tanga and Ichimarum 1974) : 


1 - <w >^/< k| 


V > 

D i 


(4b) 


where k^ is tha component of tha wave vector parallel to tha magnetic field, 
u) is tha wave frequency and the average is taken over the wave spectrum. 
Equation (4b) follows from conservation of energy and momentum of particles 
and waves. Now the dissipated power Chat goes into electrons (P^) and 
ions i?j) is 


P 

a 


5*1 "i)' 


p, "iS' 


<5) 


With for and xi * function of the drift velocity v^, the 

wave level w^ and the wave spectrum u)^(k) it is possible, in principle, to 
calculate Che above dissipation rates and P^. Because in general the 
wave level and wave spectrum are not known a priori, we must follow a special 
procedure chat is outlined in Che next section. 


2.2 Marginal stability analysis 

Recently, a method has been developed which enables one to derive values 
for Che transport coefficients of a plasma in a turbulent state, without solving 
any microscopic equations (Manhelmer and Boris, 1977; Manhetmer, 1977). The 
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oMthod us«« th« fact that an unatabla plasma avolvas toward a marginally stable 
sta"« and has already bean success fully applied to computing laboratory plasma 
pr...les (Kalfsbeek, 1978; and references cited in Manheimer and Boris, 1977). 
In this paper we shall use the marginal stability analysis only for computing 
the plasma resistivity and its essence is as follows. Suppose we start with a 
small currant so that dissipation is governed by Coulomb collisions only. 

From an analysis of the linear dispersion relation (see Appendix llA and for 
example Kindel and Kennel, 1971) one can find the minimum drift velocity for 
which the current parallel to the magnetic field becomes unstable. One sees 
(Figure lA) that for a plasma with T^/T^ < 8, ion-cyclotron waves have the 
lowest threshold, the critical drift velocity is denoted with v^^. Therefore, 
if the current density grows, the plasma will first become ion-cyclotron un- 
stable. As long as the ion-cyclotron waves are not saturated, the drift 
velocity Vj^ remains approximately equal to v^^ and thus it is possible to 
calculate the effective resistivity from 

' ^IC " (6) 

Prom this and from Tigp ■ one can calculate 


Figure 1A,B 

The energy density of the waves w^^ can now be computed from an expression 
for as a function of the wave level (see Appendix II3, formula (B7)). 

Next compare this wave level w^^ with the saturation wave level w*^ for Ion- 
cyclotron waves. If WjQ the assumption of marginal stability Is correct. 

As soon as the waves saturate, the current density starts to deviate from 
JlC for Ion-cyclotron waves and Increases till the Instability with the next 
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Iow«fC threshold mcs In. This is Chs ion*scousCic instsbillty. Again chc 
plasna avolvas coward Cha uarginally scab la scaca, this Cima with J raoaining 
elosa Co for ion-acouscic wavas (saa Figura lA). 

Apart froo cha rasiscivicy, cha narginal stability analysis also proviUaa 
cha fraquancy and Cha wava vaccor around which Cha wava spactrum is paakad and 
thus anablas ona Co calculaCa or 'lA') from aquation (4b). For 

Cha datails of this wa rafar to Appandix IIA and cha rafarancas citad Chara. 

A graph of xi < function of T^/T^ is displayad in Figura IB^ If Cha ion- 
eyclocron wavas ara saturacad, calculation of xxc ^^comas more difficult (saa 
Appendix IIC). For this case, cha calculation of vq follows from Figure 2. 

From the above discussion it follows that cha heating process consists 
ot a number of successive stages. The swicch-over from ona stage Co Che next 
is daCerminad by Che threshold and the saCuraCion level of Che waves. We 
have summarized this in a block diagram, sec Figure 2. It is unforCunace chat 
the wava levels cannot be determined very reliably, despite the decisive role 
they play in Che computation. 

Figura 2 


2.3 Ecu ac ions 

ElacCron-and ion-haating and the growth of the electric field are described 
by Che following equations: 

dT T — T 

J • - J ♦ t \ Xi nt (7a) 

•q ^ 

dT T - T 

ynkB • jnkg -h c : (l-Xi)nij^ C7b) 

•q ^ 


(7c) 
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with 



( 8 «) 

(8b) 


Is th« slsecron-lon equilibration tint (Slvuhkln, 1966)^ V 1* growth 
rata of tha alactric flald. Ua assuoa y Co ba a constant. Tha first tarm 
on tha right hand slda of tha tamparatura aquations dascrlbas hast axchanga 
batwaan alactrons and Ions by classical Couloob collisions. Nota that it doas 
not contain a factor c bacausa classical aqulllbratlon takas placa avarywhara 
Mtslda tha currant layars. 

Wa now Introduca the following dloenslonlass varlablas: 



(Uc) 
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with 


EF 


I Hi 
1 


XgF ^ I /'’EfJ Xi 


( 12 .) 


n^j, ■ (6/a) 


(m/M) 


(v^/vjj) E/T^ 


(I2b) 


a I 


T (O) 


3it 


a 

M 


^=D 


; Ep(o) 


4tt n< In A 

k,T 

B 9,0 


(12c) 


(12t) is Juic ths rtlscion n^p-E/j from (8s), rfcxprcsscd in the 
diasnsionlsss units (9). 

It so hsppsns that ths diasnsionlsss constant sr can bs vrlttsn in tscas of 
ths ratio of ths Initial alactric fiald and the initial Draictr fiald Eq(o). 
All Tl's ara axprassed in units of T1 cl(o), tha classical resistivity at t ■ o. 
Tha growth rata -y in aquation (11c) is dimensionless and expressed in units 
of t"^(o). V, is the electron thermal speed, ■ (kgl^/a)^. 

III. RESULTS 

In this section we shall present some scenarios of fast plasma heating 
computed from eq. (lla-c). At t ■ o we always have equal alectron and ion 
temperaturas . We shall start with an investigation of the classical heating 
phase. Next we will choose a specific set of parameter values and show 
numerically how the plasma evolves in the other heating stoges as well. 


3.1. Classical heating 

In stage 1, xgy " ^ (dimensionless units ). Addition 

of (11a) and (11b) gives 

~ (T,+ T^) - } E - exp (Yt) (13) 

It is easy to- solve eq. (13) for two limiting cases. The first case is strong 
electron-ion equilibration which keeps and Tj^ equal. Tne second iiaicir.g 
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cu« h«pp«ns if ch« crwiftr of <incrgy from electrons to Ions by Coulomb 
collisions is negligible. In both ceses the tempereture growth is explosive 
(Coppi end Friedlend, 1971). More interesting then the tempereture growth 
itself is the behevior of the temperature ratio during classical heating. 

We found that the explosive temperature growth does not automatically lead 
to large The reason for this is that before T«/T^ starts to deviate 

considerably from unity, the plasma becomes electrostatically unstable. The 
larger the electric field is the lower the electron and ion temperatures at 
which the plasma becomes unstable. Therefore (Tf/T^}^^, the value of T,/T^ 
at instability onset, decreases with increasing growth rote of the electric 
field. In Figure 3 we have plotted as a function of the growth rata 

of the electric field for three different filling factors c. The instability 
is always the ion-cyclotron instability. Obviously (T«/T^)cr 1* largest for 
Y * o and c * 1, when it reaches the value 3.1. Therefore it is impossible 
to make T^/T^ > 3.1 by classical Coulomb dissipation in a plasma with an 
exponentially growing electric field. As we shall show in the next section, 
the subsequent evolution of the plasma after instability onset, is not necessarily 
toward large T^/T^^. This is due to the fact that ion-cyclotron waves preferen- 
tially heat the ions. Therefor, we cannot confirm the conclusion of 
Rosner et al (1978) that a large T^/T^ can be easily produced. They reached 
this conclusion by their assumption that after the classical heating phase 
the tvo-atream (i.e. Buneman) and ion-acoustic instability will be excited. 

This is not true becajoi the ion cyclotron instability has a lover threshold. 

3.2. Results Includinz anotaalous heating (no inertial resistivity) 

We shall now present some scenarios of fast plasma heating. With this 
aim we choose the following dimensionless parameters: 
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a ■ 5 ; E • 0.01 ; y • 10 ( 14 ) 

Eqs. (lla-c) ipend at first sight only on the constants specified in (14). 

In reality, the density, magnetic field etc. all enter because the marginal 
stability analysis requires us to compute wave energy densities and to switch 
stages on the basis of this. Therefore we extend (14) with parameters 
appropriate for impulsive heating as it is assumed to occur in solar flares 
(Brown et al., 1979; Smith and Lilliequist, 1979); 

n ■ 10^^ cm“^ ; B « 500 0 ; T » T, « 5 10® K (15) 

e , o‘ 1 , o 

The density and the magnetic field are kept constant; the Coulomb logarithm 
is also kept constant, equal to 20. The above parameters imply a classical 
heating phase first, because at t ■ 0 we have * 0.05 and vj)/vj^ ■ 0.02. 

With (14) and (15) we find 

* 0.003 ; Sj^(o) « lO"^ V cm"^ 

We terminate the computation as soon as T^ “ 100 (i.e. 5 10° K) . The growth 
time of the electric field is Tg(o)/v or 0.7 s. In all scenarios considered 
below, increase of T^ by two orders of magnitude takes place in approximately 
7 or 8 growth times, that is in about 5 s, while increase of T^ from 5 10^ to 
5 10® K happens in 1 or 2 growth times, or about Is. We defer all discussion 
on application to the flare heating problem to the next section. 


Figure 4,5 
on left page 


Figure 6,7 
on opposing right 
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In stAi?« 2 th« pr«s«nc« of ion-cyclotron waves keeps the plasma in a 
marginally stable state. From Figure 13 we see that ion-cyclotron waves pre- 
ferentially heat ions. Therefore they tend to decrease the ratio while 

classical resistivity still tends to increase If the resistivity due to 

ion-cyclotron waves starts to dominate classical resistivity « the ratio 
Te/T^ will drop. Whether this will happen or not depends on the saturation 
level of the waves. Unfortunately, this saturation level is not well known. 
Therefore we consider three different saturation levels (see also Appendix IIC). 

(1) a low saturation level due to plateau formation (Petviashvili . 1964), 
strictly valid for an infinite homogeneous plasma only. In this case the 
maximum contribution of the ion-cyclotron waves to the resistivity is 

of the order of the classical resitivity. 

(2) a high saturation lev's! due to ion resonance broadening iPalmadesso et al.. 
1974). From (C2) we find that w*^ is approximately equal to 1, sometimes 
even larger than 1. For this reason we preferred to simulate saturation 

by ion resonance broadening by taking a fixed w^^ - 0.5. In this case, the 
resistivity due to icn-cyclotron wav*es soon dominates classical resistivity 

(3) An intermediate saturation level, ■ O.Cl. 

Note that all our wav'e energy densities include both the electrostatic and 
the kinetic part. Therefore, because the resistivity is determined by the 
electrostatic energy' density only, equal total energy densities in lon-cyclotror 
and ion-acoustic waves can give resistivities that are quite different. The 
wave levels for ion-acoustic waves that occur in our computations, are always 
low, so that they will not saturate. We will not pay much attention to the 
behaviour of the electron and ion temperatures themselves as a function or time. 
We found it to be much more interesting to follow the ev-oluticn of the plasma 
in Ficure lA, the drift velocity diacran. Though a point in this diagram does 
not completely specify t.he state of the plasma, it shows which wave ti-pe* 
are present and how the energy is distributed over electrons and ions. A 
part of the evolutionary paths coincides with the marginal stability curves 
(It is somewhat like the Hertcsprung-Russell diagram with evolutionary paths). 
The drift velocity diagram snows the stages, defined in Figure 2, through 
which the plasma gees during the heating process. For the three saturation levels 


mentioned above, we have calculated the evolutionary paths in the drift velooity 
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-c«t« 1- w« illustraCf the Ctoparaturts T, and , Ch# fraction \£p of tha 

dlaslpatad anargy going into alactrons, tha raalstivltlas Ipp) and 

tha wava lavals (w , ) as a functior. of tiaa (Flgura 5) . In Figura 4B tha 

plasnai;cas through staga 1 and ands in strga 2 with 0*15. 

Tha ion-cyclotron wa\-as that doainata tha heating procass apparantly never 
saturate. In Figure 4C tha plasaa goes through staga 1, axcitas ion-cyclotron 
waves (staga 2) that saturate (staga 3) and stays in staga 3 till tha and of 
tha cotaputation (T^ ^ 100); than ^ l.l. Tha reason for to remain 

close to unity m staga 3 is that tha fraction of tha dissipated energy' that 
is put into tha electrons by tha waves increases the further the drift velocity 
is above the ion-cyclotron critical drift velocity v^^ (Appendix IIC, eq. (C3) ) . 
In both Figures 4B and 4C tha increasing resistivity due to ion-cyclotron waves 
prevents the fast growth of the current density to tha ion-acoustic threshold. 

In Figure 4A ion-cyclotron waves saturate at such a low level that they are 
not important for tha final evolution of the plasma. Tha growing electric field 
drives tha plasma lon'^acoustlcally unstable (stage 4) and the plasma evolves along 
tha marginal stability curve for ion-acoustic waves to the right in tha drift 
velocity diagram. That is the electrons are preferentially heated and we end 

with T /T, \ 6.3. The fact that T /T, stabilizes at about this value when ion- 

e i e i 

acoustic waves are dominant can be easily explained (Kalfsbeek, 1978) . With 
eqs. (lla,b), ignoring all terms except the icn-acoustic one. we find 


T /T ■ X / ( I - X ) 
e'^ i ^lA ' 'lA • 


(17) 


With the help of equation (A91 and the introduction of x ■ 
for X > 6: 

V - 


dx/d(ln T,) • X - (m/X) * exp(0.5x + 1.5). 


(1C) 


It follows that X always approaches Che value for which Che r.h.s. of eq. (13) 
vanishes, chac is T«/t^ approaches che limic 6.>». This value only depend'^ on 
M/o. If energy losses are Included, chis value 6.4 will probably become lower. 
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3.3. Intrtlal restactvttv 

We aow went to discuss how the size of the current layers influences 
the heating process. We assuce that the electrons outside the current layers 
have a Maxt>rellian distribution. The electrons that enter the current layers 
with a (thermal) velocity parallel to the electric field are decelerated 
(they lose energy), while the electrons with a (thermal) speed in the opposite 
directi'tn are accelerated (they gain energy). Because the gain is somewhat 
larger than the loss, the net result is an increase in the energy of the 
electrons. The directed energy they gain in this way is randomized in between 
the current layers by Coulomb collisions and by scattering off chaotic mag- 
netic fields. Therefore, the effect of inertial resistivity is to heat the 
electrons in volume V. Of course, inertial resistivity will only be important 
if the time during which a thermal electron crosses a current layer is of 
the order of, or less than a typical (effective) collision time. This means 
that current layers with a characteristic length scale L (see Appendi.x IIB) 
smaller than the effective mean free path of an electron will carry a 
current that is determined mainly by inertial effects. Inertial resistivity 
heats only electrons (to order (m/M)^) and hence its effect is to increase 
T^/Tj^, offering the possibility of quite different paths in the drift velocity 
diagram, see Figure 6. Three different length scales will be considered, 
namely L ■ 1, 3 and IG km. For the L ■ 1 km case the influence of inertial 
resistivity is most obvious, because here T^/T^^ has increased till Just above 
8 while still in stage 1 so that ion-acoustic waves are the first to be 
generated and the saturation level of ion-cyclotron waves is Irrelevant. For 
L ■ 10 km, w^^ • 0.5, inertial resistivity is les/ dominant and ion-cyclotron 
waves are generated first. However, Instead of moving along the ion-cyc lotron 
marginal stability curve to the left as in Figure -iB without inertial res i -:tivity , 
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Che plasma now avolvas Co Cha cighc along Chis cucva. For Che inCermediaCe 
case L ■ 3 km, w|q ■ 0.01, again lon-cycloCron waves ate generated first. 
However, they quickly saturate and soon ion-acoustic waves are excited, for this 
intermediate case, inertial resistivity is also important in stage 4 as is illustrated 
by Che fact chat increases above 6.4. In Figure 7 we also show Che 

teaperacure, x£y« resistivities and wave levels for Chis case. 

IV. DISCUSSION 

4.1. Suamarv of behavior of T^/T^ during the heating 

Thou^ in principle classical heating could make » T^, chis never 
happens because the current becomes unstable before attains high values. 

The maximum value that we found at instability onset was 3.1. There appear 
to ba two possibilities for the production of a plasma ’.nth » Tj_. The 

first is that the ion-cyclotron waves saturate at a low level, so that ion- 
acoustic waves will be exci>.ad. Then T^/T^ will stabilize at 6.4. The second 
possibility is chat because of Che small dimensions of Che current layers, 
inertial resistixdty is important. In that case, large T^/T^ can be obtained 
even for high saturation levels of the ion-cyclotron waves. 

4.2. Application to solar flare problem 

If Che observed hard X-rays from the solar flare plasma are interpreted 
as thermal bremsstrahlung, it is readily calculated chat a plasma with an 

Q 

electron temperature of about 5i<i0 K (corresponding to 43 keV) and a volume 

V — 10^^ cm^ (corresponding to an emission measure 'n-dV ■ 10"*^ cm*^ when 
11 

n • 10 cm""^) is necessary. This hoc plasma must be produced fast (in a 

few seconds) because heat conduction will prevent a slow temperature rise 

fi 

to 5>1C X. Therefore, for all cases that we have presented, we used a fixed 
growth time of the electric field of 0.7 $ so that we get just about the 
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required heating titae. Not all the possibilities treated in this paper are 
equally suitable as a mechanisn for impulsive electron heating. 

T,/T^ increases considerably above unity, we expect that the neglect of ex> 
pension losses can be justified during a large part of the heating process. 
Otherwise, shock formation and disruption of the tearing mode structure 
mi|^t occur. Also, the anomalous heat conduction front theory (Brown et al, 

1979; Smith and Lilliequist, 1979) has been carried out only under the 
assumption T^/T^ » 1, so that it is uncertain if in the case T^/T^ .s. 1 the 
thermal model emits hard X-rays more efficiently than the beam target models 
(Brown, 1971). If it does not, the main advantage of the thermal model 
namely the gain in radiation efficiency, disappears. In 4.1 it is summarized 
under what circumstances a plasma with Tg » is produced. 

4.3. Exponentiating electric field 

A drawback of our work is the absence of any coupling between heating 
and the evolution of the tearing mode. One might even argue that a factor 
100 increase in T. in approximately 5 s as. illustrated in the numerical examples is 
meaningless because an exponentiating electric field eventually produces any 
desired heating. This is true but it must be realized that the largest 
electric field chat we used is only about 570.0“^ V/cm, which is very small 
compared to the maximum electric field available, The latter can be 

estimated as ■ v x B/c — 170 V/cm, where v 0.1 B/(4TrnM)^, B ■ 500 G 

—max — — 

and n ■ 10^^ cm"^. In practice will be scalier because of neglect of 

angular factors, but it will certainly be much larger than SiilO”^ V/cm. 

Apparently we have quite some freedom to use larger electric fields. We may 
need these larger fields to produc? I 5^10^K if losses are included or 
if^ for examp 1^ e is much smaller thar. 0 "“l. 
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In our calculations wa.did not taka Into account alactron runaway. Comparing 
E with tha classical Draicar field Eq (formula 12c), wa find that at t* 0 
tha ratio E/Eq «« 3)(10~^ whila at tha and of tha run E/Ep 500. Howavar, 
whathar or not alactron runaway occurs in tha currant layars dapands on an 
affactiva Draicar field vhich differs from the classical Dreicer field 

in chat the classical collision frequency is replaced by an affective colli* 

Sion frequency Vgj, ■ (ne^/m)TIgp. With Ep gy " (mv^/a) it is easy to show 

that E/Ep^gP ■ vp/v^ < 1. Neglecting electron runaway is justified 

along most of the evolutionary path in Figure 4. In the case that inertial 
resistivity is important. Figure 6, tha small size of the current 

layars prevents electrons from running away. 

4.4. Slowly changing total resistivity 

Tha behavior of the total resistivity (c.f., Figures 5B 45*2 7 b) contradicts 
the popular view in that it jumps discontinuously by many orders of 

magnitude. The essence of Che marginal stability analysis is that the time 
scale for changes in the affective resistivity (and all ocher transport 
coefficients) is reduced to the MHD evolution time. In our case, the final 
resistivity is between one and two orders of magnitude larger than at C • o. 

Of course, Che classical resistivity decreases a factor 1;(10^ since 

or T,“^*^ and T, increases by 100. Hence one might say that a factor 
10^ to 10^ in resistivity is gained. Note chat this happens only in the 
current layers, not in the bulk of the plasma. 

ACKNCWLEDGMENT 

The authors bene fitted considerably from chei.r participation to the 
CEC AM Workshop on Physics of Solar Flares, where this work has started. In 
particular we acknowledge helpful discussions with Drs. 0* S. Spicer and 

J. HeyvaerCs. 



- 19 


appbndix I 


H«at conduction around current liivrs; fmoTaCur* tiocroov and hoaogenciey . 

With th« parvMttrs that w« us« for th« tolar flare (a 10^^ cm*^; 

Tg ~ ~ 5)tl0^ la the initial phaae, tee (1S))» ve find that the mean free 

petit for eleetront at well at Iona It about 10 km. Becaute w« expect the 
current layera to have a typical trantverae length scale of a few kllooetert, 
we outt attuoe that beat tranaport from Che dltalpatloa reglont la non«elaaalcal. 
Ion heat tranaport parallel to the nagoetlc field will Chen take place with 
epproxlnataly Che Ion theroel apeed (convection) which la a hundred klloatcCera 
per second or nore > 5 x 10^). This is enough to explain the lon-teoperature . 
uniforalty parallel to the coagnetic field during the heating. Zn the aaM 
way we can argue that the electron temperature will be unlfom along Che 
■agnatic field. Perpendicular to the magnetic field the situation is more 
difficult. The classical heat conductivity la reduced considerably for both 
Che electrons and the ions. However , In Che presence of Ion-cyclotron waves, 
the Ion perpendicular heat conduction Is strongly enhanced (lonson et. al., 1979). 
Thus, within the current layers where Ion-cyclotron turbulence Is present, ion 
perpendicular heat conduction is relatively efficient, Hote that, beat flow 
awiy froQ the current layers is primarily along the magnetic field. However, 

Che magnetic field structure itself is probably quite chaotic (e.g., through the 
effects of tearing instabilities) thereby enabling the particles to distribute thermal 

energy throughout the volume V by simply following the field lines. Dif rences 
InTi^ and as well as in and T.x may in principle develop. For 
instance ion-cyclotron waves increase rather chan (Dakin, 1976, 
lonson, 1979). However, temperature anisotropies will drive fast electro- 
magnetic Instabilities (Davidson, 1972 Ch. 10,11,12; Davidson and Ogden, 1975), 
that maintain isotropy, both far electrons and ions. 
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APPENDIX XI 


A. Th« calculation of Xj^ And Xj^ undar marginally itabla conditions. 

The critical drift valocity is obtainad as follows. With the linear 
dispersion relation^ the growth rata *7 of tha waves can be datermined as a 
function of and of tha wava vactor jc. For small T is negativa for 
avary k and all wavas ara dampad. With incraasing v^, tha growth rata y will 
incraasa. Tha v^ for which tha maximum of 7(k) aguals zaro is tha critical 
drift valocity. Togathar with the corresponding wave vector , it can ba found 
from: 



If T(/T^ Is not too different from 1, an approximate expression for the 
growth rate is (Lea, 1972): 


^ “ft ^ 


Y - 


»-k|i 


(A2) 


where 


h)(k) - n 
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1 - G+ t^/t^ 


u 2 2 
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r - a I (u) > u ■ 5 
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M ' ^e “ o 
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G(u) ■ r^(y) 
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p ■ M ~ n 


/2 k,| 


is tha modified Bassal function of order n. The condition ^ • o gives 
(Laa, 1972): 
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/2 (1+ t)P 




(M) 


Miniolzlng the above expression with respect to p end ^ is equivalent to 
(Al) and in this way v^q is found. The values of A and p for which (A3) is 
fflinltulzcd can be used to calculate Xjc (^b): 

Xjj, - 1- (2m/M)‘* (T^/T^)‘*0(1+ l/A) (v^/Vj^ ) (AS) 

For A we found the following fit 

f0.t(5T /T, -J.t)'’ T /T. > 1. 

4.1 • 1 *1 

' O.Ol + 0. 16T^/T^ ^ 


The values of p chat we used, can be retraced with Che help of (AS), the fit 
for A and Figures lA and IB . 

For ion-acoustic waves we took for T^/T^ < 6, Che critical drift velocity 
v^A es given by the graph of Fried and Gould (1961). This graph is also 
reproduced in Kindel and Kennel (1971). For larger values of T^/T^ we used 
an analytical expression (Krall and Trlvelpiece 1973, Chapter 9): 


V /v - + (T exp(-0.5T^/T^ - 1.5) 

for T /T > 6 

9 X 


With the help of (4b), Kindel and Kennel (1971) and (A6) we find: 
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XA 


(l + (o/M) 



for T^/T^ ^ 1 


(A7) 




exp(-C.5T /T -1.5) 


e' 


•' ‘i 


for T /T. > 6 

e 1 


(A8) 
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wh«r« «• also uaad that ^ (n/M)^ for larga T«/Ti. For T«/T^ betwaaa 
1 and 6, an axprasslon for xia obtalnad by linaar Intarpolation. 

B. Railitlvltias 

Tha following expressions for classical and Inertial resistivity have 
been used: 

(a) For classical resistivity (Spltzer, 1962): 

'’cl" In A (Bl) 


(b) For Inertial resistivity we start from the general expression 

(Papadopoulos, 1977), where x, ■ (4tt n e^/m)^ Is the electron 
plasma frequency. Usually r Is the characteristic scattering time, but now 
we set T equal to L/v^, the average time that an electron stays Inside the 
current layer of characteristic dimension L (Spelser, 1970): 


m e 

'’in " 2 L 

ne 

Note that If the current layer Is very elongated, L Is of the order of the 
transverse dimension of the current layer. The reason Is that B will cross 
the current layer at some angle (E uid B are not necessarily parallel), 
so that a particle will cross the current layer by just following a field 
line . 

The resistivity due to electrostatic wave turbulence Is derived from 
the relation (Tange and Ichlmaru, 1974, Papadopoulos, 1977): 


(B2) 


w n 

e 3 e 




( 2 ^) 


(B3) 


• }C|| < E )c# u > Xz3 
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vh*r« altctron Ocbyt lan^th and 


“V. 




(2T) 


du < E jkf u 


(B4) 


ia tha alactroitatlc wava anargy danalty normalizad to nk^T^. 

Thaalactronic auacaptibtlity x. i* givan by (Ichimaru, 1973; Oavldaon, 1972); 


lax, ■ 


w di 2 

r 

1 

d^v 

nn 

3 

+ k. -i- 


Q V 


3v^ 

^ 5V, 


fa(v)Jn i (k|| Vj| -ui- ftfl) 

(B5) 


Tha argumant of tha Bassal function i* k^v^/D. Wa now oiaka tha following 
asaufflptlons: 

(1) only tha n ■ o rasonanca contributaa bacausa for lon-cycIotron wavaa 
kj^ Vj^ /n << 1 and for ion acoustic wavaa kj^ ^ o 

(2) f^(v), tha alactron valocity diatribution function ia a drifting Max- 
wallian with parallal to B. 

(3) |«-k,, Vjjl < k,,v^ 

Than wa find for both wava typaa: 



2 

For tha plaaoa in a marginally atabla atata wa approximata < (k,. v -u)/(k X ] > 

2 o a 

aa (k|| v^ - ui)/(kX^) , whara ky and u ara thoaa of tha marginally atabla 

wava. Tha factor (kX^)"^ ia awkward, in particular for ion-acouatic wavaa 

which hava fcrrally k ■ o for tha marginally atabla wava (Triad and Gould, 1961). 

-2 

Ha alininata (k a^) by using total (alactrostatic + aackanical) wave energy 
dansitiai. Without proof wa mantion, for k << l: 
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(kx )"^ [1 ♦ sr u - c+r,/4^)] for ic 

2(kX^) ^ for lA 

(c !• th« di«l«ctrlc eoofflciont). Thia chan laada Co: 


•r- (wc) 

0U 


(C) '’ic 



'IC 
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1 +~(i -G + r,/a^) 

^1 ‘ 
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lA 




''lA ^0 


w) 


(B7) 


(BS) 


tdiara and ara total antrgy danaitiaa. Finally, kjj - oi la 
rawrlttan In taraa of x> a* followa: 



for ZC 


for ZA 


(B9) 


(BIO) 


C. Saturation of ion-cyclotron vavta 

Zn an infinita hoooganaoua placaa, aaturation of ion-cyclotron wavaa occvira 
at a lew lavaldua to tha foraatior. of a plataau in tha alactron diatribution function. 
Tha contribution of tha aaturatad ion cyclotron wavaa to tha raaiativity than ia 
(Patviaahvili, 1964} : 
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IC 


'CL 


(Cl) 
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Bowa*/«r, duA to spAtial inhOTO 9 «n«ity of tht plaiaa, it is quits 
postibls that ion rssonancs bsoadsninq will bs tbs non-linear procsss that 
saturatss tbs instability, in that cass tbs saturation Isvsl (slsctrostatic * 
kinstie) is much highsr and is givsn by (Palmadasso at al., 1974) t 



I + rS- (1 -C+r,/d^) 

^i ^ 



(C2) 


Ms assuBS that tbs wavs Isvsl corrssponding to n._ from (Cl) and tbs rssistivity 
corrssponding to (Q2) can bs calculated with the help of equation (B7) . 

Let us assutss that in the saturated stats the wavs spsetrun is still psaksd 
about ths wavs vector where the growth rats of fonsula (A2) has a aaxisuB. 
Introduction of ths variable y ■ u/(»T ICi v ) in (A2) , so that 

, , II c 

p ■ (M/n) (A/(l +d)) (T /T, ) y , teaches us that for a given T /T. and ths 

3 si' ’ e i 

derivative ^ (y/Q) is a function of y only and doss not depend on v /v . 

r 3 - “1 “ * 

Because for the aaxiBua growth rate I (V/fl) J ■ o, the value or 


'max 


y ^ ^ does not depend on Vj^/v^, or, equivalently, (u/kjj v^) is a consta.nt for 


fixed. In the saturated state ue thus have, starting from (4b): 


where the constant has been found by requiring * Xxc * '^XC* 

Formula (C3) implies that as vq increases above vxc a progressively larger 
fraction xfc turbulently dissipated energy goes Into the electrons. 



- 94 - 


RETERENCES 

3rovm, J.C.: I9?l, Solar Physics 15 , 489. 

Brovm, J.C., Melrosa. D, and Spicar, D . S.: 1979, Ap. J . 228 , 592. 

Coppi, B. and Friadland, A. 3.: 1971, Ap. J . 1S9 , 379. 

Dakin, D.R., Tajina, T. , Banford, CJ. and Rynn, K. : 1976, J. Plasma Phys. 15 , 175 

Davidson, R.C.: 1972, Methods in Nonlinear Plasra Theory , Ch. 8,10,11 and 12, 
Academic Press (New York and London) . 

Davidson, R.C. and Ogdon, J.M. : 1975, Phys. Fluids 1045. 

Fried, B.D. and Gould, R.W.: 1961, Phys. Fluids 4_, 139. 

Ichiaaru, S.: 1973, Basic Principles of Plasaa Phvsics , Ch. 3 and 9, w.A. 
Benjamin Reading, Mass. 

Xonson, J.A., Ong, R.S.B. and Fontheim, E.G.; 1979, Planet. Space Sci. 27 , 203. 

Kalfsbeek, H.W. : 1978, Rijnhuizen Reports 78 -113 and 78- 114 . 

Kaplan, S.A. andTsytovich, V.N. : 1973, Plasma Astrophysics , Pergamon Press 
(Oxford) . 

Kindel, J.M. and Kennel, C.F.: 1971, J. Geoohys. Res . 76 , 3055. 

Krall, N.A. and Trivelpiece, A.W.; 1973, Principles of Plasaa Physics (.Mew York 
McGraw-Hill) . 

Lee, K.F.: 1972, J. Plasaa Phys. 379. 

Manheimer, W.M. and Boris, J.?.: 1977, Consents on Plasaa Phvsics and 
Controlled Fusion 3, 15. 

Manheiaer, W.M.; 1977, Phvs. Fluids 20 , 265. 

Palmadesso, P.J., Coffey, T.P., Ossakow, S.L., Papadopoulos , K. : 1974, 
Geophysical Res. Lett. 105. 

Papadopoulos, K. : 1977, Rev. Geoohys. So. Phys . 15 , 113. 

Petviashvili, V.I .2 1964, Sov. Phys. -JETP . 18 , 1C14. 

Rosner, R. , Golub, L., Coppi, B. and Vaiana, G.S.: 1978, Ap. J . 222 , 317. 

Schnack, D. and Killeen, J. : 1978, Theoretical ana Cov.putatlcnal Physics 
(IAEA, Vienna), p. 337. 

Slvukhin, D.V.: 1966, Reviews of Plasaa Physics , Consultants Bureau, 

New York. vol. 4. 

Staith, D.F. and Lilliequist, C.E.: 1979, Aa. J . 232, 532. 

Speiser, T.W.; 1970, Planet. Soace Sci . 18 . 613. 


Spicer, D.S.; 1976, NRL report 8036 . 



Spicar, D.S.: 1980, to tppaar in Solar Riytlcs . 

Spltrar, L.: 1962, Phvates of Fullv lontgad Gasa_s , Ch. 5, Intarsclence, 
York. 

Tange, T. and Ichimaru, S.: 1974, J. Phys, Soc. Jaoan . 3^, 1437. 



FIGURE CAPTIONS 


Figure lA; 


Figure IB ; 


Figure 2 ; 


Figure 3 ; 


Critical drift speeds for Ion-cyclotron waves ) and lon-acoustl< 

waves (v,-) as a function of the tenperature ratio T /T. . Also, v 
XA ® i 

v_. are nomalizecJ to the electron thenial speed v ■ (k. T /m)^ . 
lA e B e 

The nuiabers 1, 2, 3 and 4 denote the different heating stages (see 
also Figure 2) . For T /T. < 8 ion-cyclotron waves have the lower 
threshold while for T /T. > 8 ion acoustic waves are the first to go 
unstable. 

Fraction of the turbulently dissipated energy that goes into the 
electrons for ion-cyclotron waves ion-acoustic waves 

(X ) under marginally stable conditions. 

XA 

Computation of the drift velocity and the component resistivities in 
the four different plasma stages. An affirmative answer to the 
question posed in each stage, except the fourth, implies the switch- 
over to the next stage. The diagram pre-supposes < 8 and does 

not continue beyond saturation of ion acoustic waves. For T^/T^ > 8 
in stage 1 the diagram is different in that ion-acoustic waves are tt 
first to be excited. Inertial resistivity is included by replacing 

n_ by (n_. + n_„) . The symbols and expressions used are explained 
_C£s CL IN 

In the main text and Appendix II. The superscript s denotes the 
value at saturation. 

Note that for the calculation of T|'s and w's at a given time, 
the previous history of the plasma is immaterial. 

The value of the electron to ion temperature ratio when the current 
first becomes electrostatically unstable, is shown as a function 
of the growth rate of the electric field y, for three different 
values of the filling factor e. (Te/Ti)^^ is found numerically 
by solving the temperature equations in the classical phase until 
instability sets in. 


Flaura 4; 


llluscratlon of Cha avolucion of the plasma in a drift velocity 
diagram for three different ion-cyclotron saturation levels: 

(A) a low saturation level due to plateau formation (Petviashvlll , 
1964); (B) w’« * 0.5; (C) w|^ ■ 0.01. The dashed lines are the 
critical drift speeds of Figure lA. The current layers are assumed 
to be large (formally L -«), so that Inertial resistivity is not 
important ■ 0). 

Figure 5 ; The low saturation case cf Figure 4A in more detail. Time is in 
units of the growth time of the electric field, which for our 
choice of parameters is 0.7 sec. 

Figure 6 ; The evolution of the plasma in a drift velocity diagram, with 

inertial resistivity dominant j* 0). Three different length 

scales L of the current layers are considered. (A) L ■ 3 km; 

0.01; (B) L ■ 1 km, w*^ ■ 0.5; (C) L ■ 10 km, w*^ ■ 0.5. Again 
the dashed lines are the critical drift speeds of Figure lA. 

Figure 7 : The L ■ 3 km, Wj^ • 0.01 case of Figure 6A in more detail. Time 

is in units of the growth time of the electric field, that is in 


units of 0.7 sec. 



















